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1. Motivation and Introduction 

Sign Problem One of the fundamental difficulties of Monte Carlo (MC) approaches is 
known as the "sign problem" . It is encountered when the functional integrals to be evaluated 
do not have a positive definite measure. It is not related to any approximations or fundamental 
errors in the MC scheme, but it describes the situation where the statistical error can become 
very large. In general, any expectation value can be written as. 



(1.1) 



where ii and O are real-valued functions of the field variables and, in general, the measure 
need not be positive: in Lorentzian QFT this measure is complex-valued and given by 
e^'S'[(/)] 2)0^ However, if changes sign, it cannot be considered a probability density (i.e., a 
measure) . 

The standard trick to avoid this problem is to modify the measure in the following way: 
— \lJi[4>\\/ J |m[0]|- Then, one absorbs the sign of in the quantity to be measured: 

^ JO[0]sign(/.[0])A[0] ^ 2^ 
Jsign(/x[(/)])/i[(/)] 

In some cases this may work. However, the random walk guided by \lJi[4>]\ is very likely to 
predominantly sample unimportant regions in phase space. 

At this time no completely satisfactory solution to the sigh problem exists, although there 
are some promising attempts including positive projection [2, 3], fractal decomposition scheme 
[4], and Berry's phase [5] and its Stiefel manifold [6]. We believe that the techniques discussed 
in the present work lay the grounds for another possible solution path. The general success of 
the methods studied here will depend on the development of more computationally efficient 
algorithms than the ones used to illustrate the simple examples treated below. 

Smoothing out the Measure The idea behind the mollification technique is to use a 
convolution in order to smooth out and filter the measure: the highly oscillatory measure 
given by e^*^'^' T)(j) is convoluted with some suitable function (called a moUifier; see Appendix 
A for more details) and as a result only an effective contribution is left (rapid oscillations of 
the measure cancel out when integrated over the slowly varying mollifier). 

A convolution is an integral that expresses the amount of overlap of one function g as 
it is shifted over another function /. It therefore "blends" one function with another. In 
mathematics, mollifiers are smooth functions with special properties, used in distribution theory 
(generalized functions) to create a sequence of smooth functions approximating non-smooth 
functions, via a convolution. 

In an application of this idea to QFT we will replace the rapidly oscillating measure 
in Minkowski space (/ = ex.p{i S [(/)]}) by smooth functions (/^ = * ^) so that we can 
run numerical simulations that are otherwise prohibitive (see A). At this stage, the only 
requirement that is made is that: J r]e{x) dx — 1^ where the integral is performed over the 
domain of t]^. The parameter e controls the approximation, the smoothness, being made, and 
in the e — > limit, the original expressions are recovered. (The meaning of this e parameter is 
that of dilating or contracting the mollifier itself, thus controlling the range over which the 
filtering is done. Note that the calculations should be done using a non-vanishing value for e 
in order to avoid the analytical expressions which are hard to handle numerically.) 

2. Mollifying Quantum Field Theory 

The technique described above (and detailed in A) can be applied to [Lattice] QFT by 
introducing the mollified path integral 
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Z,[J] ^ J| J r?,[</> - ^] e^^I'^'-^l Vcj) . (2.1) 

This form is useful to set up numerical computations. Performing the mollification before 
performing the path integral can be very advantageous. 

Mollifying the integrand (i.e., the complex exponential of the action) smooths the highly 
oscillatory integral. The integrand changes its form, from the canonical / to the mollified J^, 
in the following fashion, 

J] = e^^I^'-^l ^ J] = (rye * /)[^; J] = \ r^e[^ - e^^I^'-^l , 

i.e., the convolution with 77^ changes variables: (j) (f. Taking [functional] derivatives of Z 
(which yield Green's functions) is just the same as taking them with respect to Z^, because 
the derivative operator commutes with the mollification. 

The [Feynman] Path Integral is constructed just as before: 



Z[J] = Jl[(/.;J]2)(/.; 
Z,[J] = Jle[^; J]Dv?. 

Note that, 

.: Z,[J] = Z[J] ; (2.2) 

where standard properties of convolutions (see [7]) have been used in the last step. 

Even though the above result is analytic, when one goes to the simulations, a small 
dependence on e shows up. (See Section 3 for a more detailed discussion on this matter.) 

2.1 Importance Sampling 

The next step consists of choosing an appropriate sampling function (for the Monte Carlo 
simulation) : 

S[^;J] 



where T)W = W[(p]T}(p^ i.e., just a change of variables in order to make the computations 
more robust: it is possible to use the integrand profile in order to speed up the calculation. 
An appropriate choice for the importance sampling function, VF, is usually given by 

wM-imi [8]. 
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A little digression is in order: If all we wanted to do was to simulate a certain [Euclidean] 
QFT on a lattice, the above reasoning would be just as valid, modulo the mollification, i.e., 
we would conclude that the proper importance sampling function was W[(f)] = The 
full-fledged formula is useless for doing computations since it implies full knowledge of the 
theory being calculated. A better choice is a simple approximation: we choose a saddle-point 
approximation in order to make the expression more manageable and still keep most of the 
characteristics of the integrand. 

Thus, the importance sampling function which is used just to better guide the 

simulation, would be nothing but the absolute value of the sum over all saddle-points of the 
theory in question. Assuming we have only one such saddle-point, the well known answer is 
given by: 



w[<p] 



,iS[M 



det{-92 + r„2 + V"[(/>o]} 



iS[4>o] 



-|tr{log(-92+m2+V"[<^o])} 



= 62 



det{-a2 + m2 + V"[,^o]}) ■ (2.3) 



Using this to find an importance sampling function for a simple scalar QFT on the lattice 
requires solving for the determinant above. In Lattice QCD, the problem of the fermion 
determinant is a very time-consuming operation that must be done for every step of the MC 
calculation. 

But this is not the only computational bottleneck. The more severe one comes from the 
so called sign problem: The exponent in ![(/)] is not bounded from below, therefore one cannot 
guarantee the ergodicity of the Markov Chain underlying the Monte Carlo draws, i.e., a much 
bigger number of draws would be needed in order to yield any meaningful answer. This is 
the reason why calculations are done in Euclidean space rather than Lorentzian/Minkowski 
space. In Euclidean space (Wick-rotating the integrand), the exponent in question becomes 
— e"*^'^': it is bounded from below and the Monte Carlo method works fine (the Markov 
Chain behind it becomes ergodic). 

It is to tackle this problem that the method of mollifiers comes in: mollifying ![(/)] will 
yield a smooth function, /e[0], whose properties help the convergence of the MC computation. 
Furthermore, if this can really be accomplished and calculations with imaginary exponents 
become a reality, the next logical step is to analyze the different phases of the given QFT. As 
shown in [9] , the different phases of a theory can be picked out with a mere choice of boundary 
conditions, which is the same as properly defining the measure of the [Feynman] Path Integral. 
The upshot is that if we want to obtain answers other than the ones that can be reached via 
perturbation theory, either the measure of the path integral or the boundary conditions of the 
Schwinger-Dyson equations have to be modified [9] . 
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Once there is nothing preventing the measure of the [Feynman] Path Integral, from being 
(in the most general case) complex-valued — in order to account for the phase structure of the 
theory — , the fact that its integrand is also complex- valued has to be taken more seriously, 
afterall this is a non-perturbative result: a perturbative series only works if we know, a priori, 
in which phase we are working, so it can be tailored to that particular sector of the theory. 

Thus, the need to address the sign problem becomes even clearer: it is not simply a 
problem of doing MC in Lorentzian/Minkowski space, it is mostly a problem of being able to 
compute all possible solutions — all the different phases — of a given QFT. 

2.2 Mollifying the Importance Sampling Function 

As outlined before, the idea is to use the mollification technique to handle the highly oscillatory 
terms present in the importance sampling function. 

There are three possible choices for the sampling function W: 

1. Taylor-expand and mollify the complex integrand: no knowledge about the saddle-points 
is necessary, but it is very computer intensive, which is the reason we will not focus on 



2. As done in (2.5), where the integrand is saddle-point expanded and mollified, but the 
explicit form of the mollifier is used in order to perform the remaining integral: not so 
heavy on the computer, but some knowledge about the saddle-points is needed; 

3. As done in (2.6), where the integrand is saddle-point expanded and mollified, generalizing 
the standard saddle-point approximation: less intensive of all 3, but even more knowledge 
about the saddle-points is needed. 

The difference between the second and third methods above is that the integral in (2.4) 
is explicitly carried out in the former, but not in the latter. As explained below, the third 
method (listed above) consists of the simple generalization of the well known saddle-point 
approximation; however, we need to note that more knowledge about the particular saddle- 
points is needed in order to perform the contour integration involved in this scheme (and this 
is the difference with respect to the second method). 

The calculation below shows the second method outlined above, known as 2'^'^-order 
gradient approximation: 



Now, a saddle-point expansion is performed on the action, 
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A mollifier needs to be chosen, and for the purposes of this calculation, a Gaussian one is our 
choice: 



ri,[cf)-(p] =e" 



WA4>] = 



WA4>] = 



2 
OO 



V2 7re2 



-<^)2/e2 i (S[vo]+| (f-vo? S"bo]) 



exp 



{i5[^o] 



+ 



1 (</i-yo) 5"[yo] ((?!'-yo) 

2 l-ieS'"[<j9o]e 



det{l -ie5"[^o]e} 



(2.4) 
(2.5) 



The important point to note about the derivation above is that the explicit form of the 
molhfier, 77^, had to be used and, because of that and the nature of the MC simulation, we do 
not necessarily need a good knowledge about Lp^. This means that we can compute the above 
importance sampling function, (2.5), and use a trial Lp^^^^^\ the [MC] simulation will do the 
job of moving towards the exact (/pq and pick out the different phases of the theory. As shown 
elsewhere, [11, 12, 13, 8], this works quite fine. 

However, for more relevant cases (e.g. 4-dimensional QFTs) the time required for simula- 
tions becomes unrealistic. The bottom-line is that a simple-minded implementation of the 
method will have the code/computer doing the work of mapping the (0, (/9)-space, while one 
could use that information beforehand in order to speed things. 

The way out of this is to generalize the following well known result: J(s) = g{z) e^^^^^ dz - 
V2^g{zo)e'f^''^e'^/^\sf"{zo)\'; 5 £ R, 2 e C; i.e., 



TT 



5[vo] 



det{5"[H}| 



(2.6) 



At this point, the similarity between (2.3) and (2.6) is clear, in fact, they are the same 
formula, except that, in the former case there is no mollification while in the latter there is 
mollification. 



3. Tuning the Mollification 

At this point we are almost done with the analysis of the mollification technique, the only 
aspect remaining being the optimal choice of e, the parameter that regulates the approximation. 

From the discussion in appendix A, it is clear that in the e limit we recover the 
original theory. However, numerically, this is not exactly what happens, because of numerical 
and statistical fluctuations. Therefore, the choice of an optimal value for e is crucial. 

Using an information theoretic viewpoint, we can choose e to optimally compress We[4>\ 
around each stationary phase point. Therefore, we can define an information entropy based 
on We[0] in the following way: 
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The optimal e is the one that minimizes S[e]. This value corresponds to maximal compression 
of the information in the importance sampling function. 

In order to illustrate the above, let us consider a very simple [0-dimensional] example, 
given by S'[0] = 0^/2. In this case, the saddle-point is 00 = and the importance sampling 
function and the information theoretic entropy are given by: 



-0V2e2 iS[4>o] 



^/S^ ' 



§[e] 



me v27r 
log(me \/2 7r ) — 
1 — m log(m) 



Plotted below are the graphs of the importance sampling function (also showing how it 
varies with e) and of the information theoretic entropy (both using m = 1): 



IrifiDrtB-EwSEiinplk-ignjncUDn: 




Figure 1: Importance Sampling function (leftmost plot, m = 1, We[0] = e / e V2 tt ') and information 

theoretic entropy (rightmost plot, m = 1, S[e] = log (e a/2 tt ) — e^) for the action S[(l)] = rri^ 0^/2. 

We can clearly see the importance function peaking around = 0, which is the expected 
behavior (i.e., it happens around the saddle-point 00 = 0), while the information theoretic 
entropy shows that the best parameter for the mollification is e ^ 0: this is not desirable, 
once lim^^o S(^) — oc, but it simply shows that we should use a small non- vanishing value 
for e, as mentioned before in Section 1. 
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4. Simulation Details 



In order to perform the simulation as outlined in 4.2, some technical details had to be taken 
into more serious account. Among them, the important ones are: random-number generation 
and non-local algorithms. Those will be the bottle-necks, as noted in [14]. 

Random Number Generators Even though some consider the random-number question 
a solved one, this is not always the case [14]. For this reason, the choice of which random- 
number generator (RNG) to use is still a critical one, not only for speed reasons but because of 
systematic errors as well. In order to address this issue the choice made was for the Mersenne 
Twister RNG [15] (period of 2^^^^^ - 1). 

Non-local Algorithms As for non-local algorithms, the choice made was for a Genetic Algo- 
rithm (GA) [16]. The "population", in the Lattice QFT sense, consists of the lattice {i.e. space- 
time) points. For a given "recombination rate", an initial population gets randomly divided into 
pairs of [lattice] points and "crossing-over" (or "genetic") operators are applied to these pairs 
in order to generate different sets of pairs. In this context, the crossing-over operators exchange 
the coordinates of given points, thus there are 4 of those operators, namely: T^, T^, Ty ^Tz. 
That is, if the crossing-over operator Tx is applied to the points (j) = {(j)t^(j)x^(l)y^<l)z) and 
^' {^'t^'x^'y^'z) the outcome wiU be: r^[0,(/)'] {{(t)t,(t)'xAyAz)^ (0^,0^,0^,0'^)}. Anal- 
ogous definitions are valid for the other crossing-over operators. The key features for this 
choice were: which Ti to use is random; whether or not the change (j) ^ (j)' is made depends 
on a given probability and the recombination rate can be arbitrarily chosen (although keeping 
it below 1% showed to be a good tune). Moreover, the genetic operators are unitary, i.e., 
Ti = . This implies that = 1, which guarantees the so-called detailed balance of the 
MC simulation (this is the ergodicity of the algorithm). At this point, the only step remaining 
is the explanation of the probability profile used in algorithm 2 below. 

1: Choose the recombination rate such that: ^ r ^ 1 {usually between 0.5% and 1.0%} 
2: Draw a random number f3 G [0, 1] and compare with r 
3: if /? > r then 

4: Metropolis Monte Carlo loop {using the random walk technique} 
5: else 

6: Genetic loop {see algorithm 2 below} 
7: end if 

Algorithm 1: Genetic Algorithm 

In the notation used below, P is such probability profile and P{(j)[,(j)j) means that this 
profile is calculated using the [GA generated] points 0^ and 0^-. It proved useful to chose the 
mollified importance sampling function as this profile P. Also, note that 0^ and 0j, just like 
(/)[ and 0^-, are just particular values of the field at the sites i and j. A comparison between 
these two profiles is then performed. 
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4: 



3: 



2: 



1 



Choose a random pair: {^i, ^2} = {(/)(ti, xi, yi, zi), 0(^2, X2, ^2, ^2)} 
Generate the pair = {^(ti, x^, 2/^,2;;), x'^, y'^, 4)} = Ti{(j)i, ^2) 

Draw a random number c G [0, 1] 
if c>P((/);,0'2)/^(0i.02) then 



6: 



5; 



do nothing 
else 



8: 



7: 



perform the exchange (0i,02) ^ (j)^) 
end if 



Algorithm 2: Genetic Loop 



4.1 Euclidean Lattice Simulations 

In order to be able to understand better what will happen in the Lorentzian/Minkowski case, 
a quick and dirty reminder of the Euclidean one is presented below. 

In what follows, $ is the set of all possible field configurations and 1 1$| | is its cardinality (i.e., 
the number of elements in the set). For each element of $ (i.e., for each field configuration) one 
has to average the observable over all the lattice points (which are computed via Metropolis). 
Each element in the Markov Chain (generated by the Metropolis MC algorithm) is denoted 
by 0^. To obtain the final result, an average over all field configurations is made. (To make 
clear that the final goal is the actual lattice computation, the measure is denoted as [dcj)].) 



4.2 Mollified Minkowski Lattice Simulations 

In the Minkowski version of the above, the Wick rotation is not performed, therefore we are 
left with the original form of the exponent. 



The above functional will be mollified in order to yield more tractable expressions. The 
fact that there will be two functional integrations (rather than just one, like above) should 
not bring many problems. The real question here stems from the fact that the mollifier, 77, 
mixes the two of them. 



JO[</>]exp{-5[</>]}[#] 
Jexp{-5M}[#] 




fO[(/.] exp{i5[</)]}[#] 
Jexp{i5[(/.]}[#] 
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S{SVe[cf>-V>]0[^] exp{iS[^]} 

ll^ll ll^ll h I E&i - ^W] exp{i5Mq}/Wf[^W] J ' 

where it is understood that the configurations in $ are chosen with respect to the importance 
samphng function, VF^^^i^], and the configurations in F are chosen with respect to a uniform 
distribution. This happens because the mohification process (namely the convolution) mixes 
the variables from $ and F together. Therefore, the only way to implement this "mixing" is 
by having a uniform distribution for F and implementing the "interaction" with the variables 
in $ — via the mollification — explicitly, using 77. Note that ||$|| = ||i^||- 



(0) = 



5. Spontaneous Symmetry Breaking 

Note that in this paper we use the designation "spontaneous symmetry breaking" somewhat 
loosely in that we apply this term to all solutions which do not approach those of coupling 
constant perturbation theory in the limit where the coupling vanishes. We have been discussing 
lattice QFT in Lorentzian/Minkowski spacetime because one of the main objectives of this 
paper is to examine a numerical approach to directly calculate Green's functions in different 
phases in QFT. This necessitates being able to evaluate path integrals of complex exponentials. 
The arguments for this have been given elsewhere [9] and will be illustrated in a particular 
example in the following section, but it is fairly easy to understand why this is the case. 
The traditional [Feynman] Path Integral formulation of a QFT involves integration of the 
exponential of the action over every field variable at every space time point. Usually, these 
integrations (assume the action is written in terms of self adjoint fields) range from negative to 
positive infinity along the real axis. In the limit of small couplings, this form of the [Feynman] 
Path Integral generates perturbation theory and the results appear (for finite number of 
expansion terms) to be regular at vanishing coupling. Thus, all expansions that are not 
regular at vanishing coupling, such as the traditional symmetry breaking expansion of quartic 
scalar field couplings, are excluded. In order to avoid this restriction and produce all possible 
solutions of the QFT it is necessary to extend the [Feynman] Path Integral integrations to 
complex values of the fields in a way consistent with the field equations and reality properties 
of the theory. 

Indeed, as shown in [9] and section 6, the different phases of the theory emerge through 
the varied boundary conditions consistent with the equations of motion. Equivalently, rather 
than varying the boundary conditions of equations of motion, the measure of the Path 
Integral can be changed. In general the number of choice of paths of integration for the path 
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integral correspond to the number of independent solutions of the Schwinger-Dyson differential 
equations. 

If we define a QFT via its [Feynman] Path Integral, all we need to know is the action, 
S[(l)^ J] (for a field (j) whose source is J), of a given model, for then we can write the generating 
functional as: 



Z[J] =K Jexp{z5[(/); 



where N is a normalization constant such that Z[J — {)] — 1. The crucial question that remains 
unanswered in this approach is: ''How does one properly define the measure Vcj)?^^ 

The best answer so far (for 4-dimensional systems) says that this can only be done for 
free QFTs, via the use of cylindrical functions [17, 18]. (Note that this only happens in the 
continuum. In its Lattice formulation, QFT is free from such peculiarities because of the 
lattice regularization.) 

Analogously, defining a QFT via its Schwinger-Dyson equation, 

^^^tlilz[J] - J(x) Z[J] = ; 

is equivalent to substituting (f) —i jj in the action, requiring the equations of motion to 
be the solutions that extremize it. Note that the Schwinger-Dyson equations are a system of 
infinitely many [partial] differential equations, one per each point of spacetime. 

Therefore, when thinking in terms of differential equations, the boundary conditions are 
responsible for the phase structure of the theory. Intuitively, the picture that comes to mind is 
that of a portion of space divided into as many subsets as there are solutions to our equations 
of motion, such that in each of those regions, the equations of motion satisfy appropriate 
boundary conditions. 

On the other hand, when thinking about the [Feynman] Path Integral, usually it does 
not seem bothersome that the measure is not properly well-defined. In fact, determining 
the boundary conditions for the Schwinger-Dyson equations is analogous to determining 
the measure for the Path Integral. Thus, just like the boundary conditions, the measure is 
responsible for the phase structure of the theory (in the integral representation of the problem). 

For completeness sakes, this is how a QFT in different spacetime dimensions, behaves: 

d = For 0-dimensional QFTs — i.e., QFT on a point (the [Feynman] Path Integral degen- 
erates into a simple integral) — there is no such thing as phase transition since there is 
no such thing as dynamics. However, the phase structure of the theory survives, given 
by the different boundary conditions (resp. measure) needed in order to determine all 
the solutions to the equations of motion. 

d = 1 For 1-dimensional QFTs — i.e.. Quantum Mechanics — , again, there is no such thing 
as phase transitions, for if the theory has 2 different vacua we could take a linear 



- 11 - 



combination of them to be the "real" vacuum state, given that they would be related by 
tunneling. (For the analogous case in Condensed Matter physics, please refer to [19].) 

d ^ 2 In this case, phase structure and transition exist; this is the complete scenario. (Note 
that there is no operator that relates 2 inequivalent 0- vacua of the theory, therefore they 
belong to different algebras.) 

The real lesson to be learned from all of this is that Spontaneous Symmetry Breaking is a 
phenomena generated by boundary conditions; whether you use them to define the measure 
of the Path Integral or to define the Schwinger-Dyson equation is just a matter of personal 
preference. It is not said anywhere that the limits of the Path Integral have to be real; what 
we need to have are real observables. In fact, as shown below on section 6, it turns out that 
in order to have symmetry breaking we need a measure that is not necessarily real: this will 
enable the computation of all the solutions of a given QFT [9]. 



6. Results 



In what follows, the results obtained thus far are presented. In short, they are in 0, 1 (time) 
and 4 spacetime dimensions. 

6.1 Lower Dimensional 

Let us start by addressing the 0- and 1-dimensional results. Plainly and simply put, this 
means that we are solving a simple integral in spacetime dimensions since the [Feynman] 
Path Integral degenerates into a standard integral and, in the 1-dimensional (time) case, one 
will be doing Quantum Mechanics. 

For the 1-dimensional results, refer to [8]. There, the Quantum Chemistry of the problem 
is fully treated and addressed. Note, however, that the importance sampling function chosen 
in [8] is different than the one used in this work: the quantum chemistry was done using (2.5) 
while we use (2.6). 

Below, the 0-dimensional results are summarized and, in order to illustrate the features of 
the mollifier technique, different properties of those models are made explicit. 

Airy Function The action for this model is given by: S{x) = + Jx. This is 

an interesting model because we can explicitly calculate the partition function and com- 
pare it with the results coming from the Mollified Monte Carlo procedure. Moreover, this 
model has 2 stationary phase points and the results displayed are from the one in the com- 
plex plane which is not accessible with normal Monte Carlo. The quantity of interest is: 



Figure 2: Plot of the Airy potential, V{x) = 
highlighting the fact that it is not bounded from below. 
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X!^ exp{i^ + iJx}dx ^ Ai(J) 

The [first two] graphs below show the highly oscillatory behavior of the integrand in the 
partition function: the left one is its the real part, while the right one is its the imaginary 
part. 




Figure 3: Real and Imaginary parts of e*^ /3+z gj^Q^jj^g highly oscillatory nature of the problem. 

Below we see the graphs of a particular 2-dimensional slice of the above pair, where 
J = — 1: the real (left) and imaginary (right) part of the mollified Airy-integrand (in red, 
e = 0.1) is in contrast to the non-mollified integrand (green). 





Figure 4: Superimposed plots of mollified (red) and non-mollified (green) real and imaginary parts of the 
integrand e'^^ /3-ix ^ showing the smoothness achieved on the J = — 1 two-dimensional slice of the previous 

graphs. 



Below we see the graphs of the real (left) and imaginary (right) part of the mollified 
Airy-integrand for different values of the field x and of the mollification parameter e: analogous 
to Figure 3 but for the mollified version of the integrand (which corresponds to the red lined 
plots of Figure 4) . 
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UDlllbd IntngrarKl {Imaglnaiy Pjit, J-^i;i 




Figure 5: Plots of the real and imaginary parts of the mollified integrand for various values of x and e. 



The 2 stationary phase points are given by xq G {=b\/— J }, where xq is the solution of 
^^'(^o) = 0, and the [properly mollified] importance sampling functions are: 



»'/>°M=exp|-i?^}+exp|-l?^|. j>0; (6.1) 
=exp|-l(i^)1+exp|-l(^i±^)1; J<0; (6.2) 



It is customary to see the source term above called t and identified with time but this is not 
the approach taken here. 

Below we see the graphs of We with respect to the e parameter and x, for fixed values of 
the source: on the left J = 4, and on the right J = —16, 



ImpcnannSiSipInq Flncrtan: Ary Acllan, J-+ 




Figure 6: Plots of (6.1), for J = 4, and (6.2), for J = -16. 
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Figure 7: Plots of 2-dimensional slices of (6.1) and (6.2) for various values of e. 

It is not difficult to see that there is an optimal value for the parameter 6, as shown below 
in the graph for the information theoretic entropy (see Appendix B for further details): we 
have the information theoretic entropy for J > on the left, and J < on the right, 
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Figure 8: Information theoretical entropy for finding the optimal value of the mollification parameter: on the 

left we have J ^ 0, and on the right we have J < 0. 

Furthermore, it is easy to perform the analytical calculations outlined above (using the 
Airy functions on the Generating Functional, rather than its mollification) to get the pure 
evaluation of this model. The results show the mollified answers (given a proper choice of 
e) are sensationally accurate and show, already in this simple example, that a region not 
allowable in ordinary Monte Carlo approaches can be computed. In fact, due to the simplicity 
of this example, modern computers can calculate it straightforwardly, without the need of 
special tricks to handle the highly oscillatory integrand. The point of this example is twofold: 
show the precision and the accuracy that the mollifier method can achieve and also to illustrate 
how this smoothing procedure works. 
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0-dimensional (p^ Theory Also known as ultra-local 0^, its action is given by: S[(l)] = 
IJi(j)^ /2 + ^0^/4. As before, the saddle-point 00 is such that ^^[^o] = 0. The graphs for the 
Action above (positive mass on the left, negative mass on the right) are given by: 




Figure 9: 0^ potential for positive and negative values of /i. 

Below, we show the highly oscillatory behavior of the integrand of the Partition Function 
given by the above action, i.e., / = e^^^^'^ for positive and negative values of the fi 
parameter: 



rt(V»*^ -m^t'.T + B-j'rt - Jtt 11.1 .^1 





Figure 10: Plots of the Real and Imaginary parts of the integrand for various values of J and fi — 1^ g — 1 

and /i = —1, g = 1. 

Below, we see the real (left) and imaginary (right) parts of the mollified (/)^-integrand 
(red) in comparison with their non-mollified counterparts (green): 
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Figure 11: Superimposed mollified (red) and non-mollified (green) real and imaginary parts of the integrand 

for/i=l,5f = l, J=— 1, e = 0.1. 

Below, we see the real (left) and imaginary (right) parts of the mollified (/)^-integrand for 
different values of the field (f) and of the mollification parameter e: 
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Figure 12: Plots of the real and imaginary parts of the mollified integrand for various values of the field and 

the parameter e. 

As explained in [9, 12], the three different solutions/phases of this model can be selected 
via an appropriate choice of boundary conditions: = R, r+ = (— oc, 0) U (0, i oc) and = 
(— oc, 0) U (0, —i oc). The leftmost figure below shows those boundaries in the Argand Plane: 
= R is the green line, r+ = (— oc, 0) U (0, i oc) is the red line while = (— oc, 0) U (0, —i oc) 
is the blue line. The rightmost figure below shows the regions (shaded) of the complex (/)-plane, 
{(j) — pe^^)^ defined by cos(4 0) ^ 0. Any contour, starting and ending at infinity, within one 
of these four domains corresponds to a particular solution of the 0-dimensional 0^. 
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Figure 13: Left: Three different contours that render the Path Integral finite and depict the three distinct 
solutions obtained in this model. Right: Any contour that starts and finishes inside the same-color shaded 
regions (following the general directions shown on the left) corresponds to a particular finite solution. 

Using the three different contours above we are able to find the three different solutions to 
the equation of motion of this 0-dimensional QFT. The Schwinger-Dyson equation is given by, 

where we can easily see (again) that there has to be three solutions. These are called Parabolic 
Cylinder Functions and are denoted U{fi/g^ J), V{fi/g^ J), J). Those behave in the 

following manner (see Appendix C): 

Regular at g ^ Consistent with perturbation theory, 
Singular y^) at ^ ^ Symmetry Breaking, 
Singular exp{/x/4^}) at ^ ^ Instanton. 
Thus, 

X!^ Jro,± r]e((l) - (p) (D[(p] explil^if'^ -{-i^if"^ -iJ(f\d(fd(f) 

{0)e[J] = 4 ^ . 

iZo Jro>± ^^(^ ~ ^) ^^Pj^ 2 + ^ 4 ^^^^ 

It should be noted, though, that depending on the chosen contour (F^, r+, F"), the above 
integral will resemble a Fresnel one and, as such, its computation is quite simplified. 
For the graphs below, we used the following values for the saddle-points, 

00 G {0,±V-m/2^ } , 
where ^^f^o] = 0. Therefore, the different importance sampling functions are: 
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Symmetric Phase Given by the contour, this is the one which is regular in the hmit 
^ ^ and accessible via perturbation theory, 

l^f[0] = exp|-l(^)'|; 00 = 0; 

Solitonic Phase This phase is given by the linear combination of the contours r+ and 
such that /X ^ 0, i.e., this represents the soliton solution, 

,rf>>l=exp{4(:^^±M)!)}+exp{-l(*^^)}; = ±,, ; (m > 0) ; 

Broken-Symmetric Phase This one is given by the linear combination of the contours 
r+ and r~ such that /i < 0, i.e., this represents the solution usually referred to as 
broken-symmetric, 

lFrVl=exp{-i(*±*)'}+exp{-l(t*)'}; 4>t = ±V^I^ : (p<^)- 

Below, we have the graphs of these in imp«i«i«ainpinBFur«dDn;^*Ai:im^t^i.viirE?j,F-i.B-i 




Figure 14: Plots of the mollified importance sampling functions for various values of the field and e: 
symmetric (00 — 0), solitonic [fi — 1^ g — 1) and broken-symmetric (/i = —1, g — 1). 

The graphs below show the entropy for the above importance sampling functions, as 
discussed in Appendix B: 
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Figure 15: Information theoretical entropy for finding the optimal e in the symmetric, solitonic and 

broken-symmetric cases. 

6.2 Higher Dimensional 

As for the 4-dimensional resuhs [for a Free Scalar QFT] the news is not encouraging. As 
mentioned in section 2.2, there are severe efficiency constraints in this implementation of the 
mollification; in section 2.2 three possible choices of solutions to this technical issue were given. 

It turns out that, from those three solutions, one is non-efficient while the other two 
produce analogous results: only the imaginary parts of physically interesting quantities 
(generating functional and Green's functions) can be tamed. The real parts continue to show 
an oscillatory pattern. 

Although this is at least disappointing, the solutions obtained are good enough to control 
the imaginary part's oscillations. 

It is our hope that changes in algorithms will allow some version of mollification to work 
in real world problems. 
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A. Convolutions and Smoothing 

With the aid of the tools developed below, the smooth approximations of functions can be 
done in a mathematically rigorous fashion. This is useful to justify the statements made in 
this work, where this technology is applied to the generating functional of an arbitrary QFT. 
For proofs of the theorems shown below, see [20] [21]. 

A.l MoUifiers 

Notation. If U CMJ^ is open, dU is its boundary and e > 0, let 11^ = {x e U \ dist(x, dU) > 
e}. Further, let B{0^e) be the ball centered on and with radius e. 
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Definition A.l. (Standard mollifier.) 

A mollifier^ 77^ also called an approximate identity^ is a positive "^^(R^) function. The 
standard mollifier is defined in the following way: 



• Define r] G ^^(R^) to be, 

C exp(^j^^y if \x\ < 1; 
0, if \x\ ^ 1. 



r]{x) = 



The constant C > selected so that J^^ r]{x) dx — 1. 
• > 0, set 

77 is called the standard mollifier. The functions 77^ G satisfy J^nr]s{x)dx = 1 and 
supp(%) C B{0,£). 

Definition A. 2. (Mollification.) If f : U ^ is locally integrable, define its mollification to 
be, 

That is, y X ^ U^, 

/ Ve{x-y)f{y)dy= / r]^{y) f{x - y) dy . (A.l) 

Ju Jb{o,£) 

A. 2 Properties of MoUifiers 
Theorem A. 3. (Properties of mollifier s.) 

1. r e -^^^m; 

~^ /; cilf^ost everywhere, as e ^ 0; 
3. If f ^ ^{U), then f^ ^ f uniformly on compact subsets ofU; & 
I If 1 ^p< 00 and fe Ll,{U), then f^fm Ll,{U). 

B. Entropy Calculations 

In the cases considered in 6.1, the information theoretic entropy can be computed analytically. 
There are 2 types of functions to be considered: 



fl{x) = exp I 
= exp 



1 / I \ 2 
1 I X ^ c 



2 V e 
1 f x^ + c 
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where c is a constant scalar. The first type above deals with the Airy Function action for 
J < (c = ±V— J) and the c/)^ action with ^ < (c = ±y^— ^/2^'), while the second type 
above handles the Airy Function action for J ^ (c = zbJ) and the 0^ action for /x ^ 

The information theoretic entropy has the general form, 



Jfeix) \0g{flix))dx 

f fi{x) dx + log ( / fi{x) dx 



Therefore, we see that the analytical results are given by: 

(e/2) ^/2^ 



§2(6, C) 



\/27r e + \og{^/2^^ e) 

e--^/^^' V^(e + c/2e) 



And, thus, we have closed forms for the entropies of the examples given. The graphs below 
show these entropies: the leftmost one shows S"^(€) while the rightmost one shows S^(€,c) for 
a range of c values. 
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C. Parabolic Cylinder Functions 

The differential equation for the Parabolic Cylinder functions is usually written as, 

-^ + {az^ + bz + c)f = 0; 

where z G C. 

The 0-dimensional 0^ action is given by S[(f)] = 0^/2+^ 0^/4, which yields two Schwinger- 
Dyson equations: 
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Perturbative solution: — i 6jZ — J 2, = ; 
Non-perturbative solutions: SjZ + JZ — /3Z = 0; 

where f3 = fi/g. The match to the previous form is obtained when a = 0, 6 = +1, c = —(3. 
Note that there are two non-perturbative solutions, one being solitonic and the other one 
being called simply broken- symmetric. 
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